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Abstract 

We model the effects of confinement on the stable self-organized patterns obtained in the non- 
equilibrium steady states of mixtures of molecular motors and microtubules. In experiments 
[Nedelec et al. Nature, 389, 305 (1997); Surrey et al, Science, 292, 1167 (2001)] performed 
in a quasi-two-dimensional confined geometry, microtubules are oriented by complexes of motor 
proteins. This interaction yields a variety of patterns, including arangements of asters, vortices and 
disordered configurations. We model this system via a two-dimensional vector field describing the 
local coarse-grained microtubule orientation and two scalar density fields associated to molecular 
motors. These scalar fields describe motors which either attach to and move along microtubules 
or diffuse freely within the solvent. Transitions between single aster, spiral and vortex states are 
obtained as a consequence of confinement, as parameters in our model are varied. We also obtain 
other novel states, including "outward asters", distorted vortices and lattices of asters on con- 
fined systems. We calculate the steady state distribution of bound and free motors in aster and 
vortex configurations of microtubules and compare these to our simulation results, providing qual- 
itative arguments for the stability of different patterns in various regimes of parameter space. We 
also study the role of crowding or "saturation" effects on the density profiles of motors in asters, 
discussing the role of such effects in stabilizing single asters. 
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I. INTRODUCTION 



The mitotic spindle in a dividing eukaryotic cell is comprised of several millions of inter- 
acting protein molecules^]. Remarkably, these molecular constituents self-organize to yield 
patterns at the scale of microns. The existence of such self-organized non-equilibrium struc- 
tures at the sub-cellular scale is a common feature of biological systems. Such structures in- 
clude the endoplasmic reticulum and the Golgi complex, membrane-bound organelles which 
participate in intra-cellular trafficking. They also include the cytoskeleton, a cell-spanning 
network of polymers such as microtubules, actin filaments and intermediate filaments^]. 

he mitotic spindle consists of two interpenetrating arrays of microtubules!^. Each such 
array originates from a microtubule organizing centre formed at the two ends of the dividing 
cell. Microtubules are semi-fiexible polymers formed by polymerizing an asymmetric dimeric 
unit0|. An individual microtubule is a polar object: microtubule ends, labelled as — and +, 
grow and shrink at different rates. This polarity dictates the direction of motion of a class 
of molecular motor proteins on microtubules. Motor proteins such as kinesins use energy 
derived from adenosine tri-phosphate (ATP) hydrolysis to exert forces and to translocate 
along microtubulesyi. The direction of motion is dictated by the filament polarity - both 
plus-end and minus-end directed motors exist 

The principal microtubule organizing centre in most animal cells is the centrosome. In a 
dividing cell, the centrosome first duplicates to form the two poles of a mitotic spindle. Each 
centrosome nucleates a radial array of microtubules called an aster. The two asters move 
apart to opposite ends of the cell as mitosis proceeds. Microtubules in the region between 
them then preferentially elongate, forming the spindle [s]. 

Interestingly, experiments on centrosome-free fragments of the cytosol containing both 
motors and microtubules also obtain self-organized aster-like structures^l. Asters, in addi- 
tion to other complex patterns such as vortices and disordered aster-vortex mixtures, are 
also seen in vitro, in experiments on mixtures of molecular motors and microtubules Q]. 
Features of spindle formation are reproduced in mixtures of motors, microtubules and gold 
beads coated with DNA ^ . The fact that such experiments are able to mimic the complex 
self-organized states seen in living cells indicates that simple mesoscale models which work 
with fewer components may be useful in capturing some aspects of cellular pattern formation 
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This paper presents a theory of pattern formation in mixtures of molecular motors and 
microtubules in a confined geometry. We motivate hydrodynamic equations of motion for 
a coarse-grained field representing the local orientation of microtubules as well as for local 
motor density fields [l^]. In our theoretical description, as well as in the experiments, rnipro- 
tubules are oriented by complexes of bound motors, yielding patterns at large scales 
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The experimental work of relevance to this paper is the following: Nedelec and col- 
laborators study pattern formation in mixtures of complexes of conventional kinesins with 
microtubules in a confined quasi- two-dimensional geometry Q]. The later experiments of 
Surrey et al. investigate pattern formation in larger systems where the effects of boundaries 



appear negligible 



13|. The experiments are supple mented by simulations which reproduce 



many features of the experiments |7|,ll3|,ll5|,ll6|,ll7|. 

Our approach is closest in spirit to that of Lee and Kardar(LK) The LK model 

is a hydrodynamic description of two coupled fields. One of these is a two-dimensional 
vector field describing the coarse-grained orientation of the microtubule. The other is a 
conserved scalar field representing the local motor density. The LK model captures two 
prominent features of the experiments of Nedelec et al.: the presence of stable vortices 
formed by microtubules and the instability of an aster formed in the early stages of pattern 
formation to a single, stable vortex at large motor densities. However, despite this success 
of the LK model, several features of the experiments are incompletely understood. The 
transition between a single aster and a single vortex seen in the experiments appears to be 
driven primarily by confinement. (When a growing microtubule hits the wall of the confining 
region, it bends and distorts. The elastic energy cost of this distortion can be relaxed in part 
if the microtubules assume a vortex configuration.) In the LK model, confinement does not 
appear to play a vital role, and the single aster to single vortex transformation is a generic 
feature even for large system sizes, provided the motor density is sufficiently large. 

Experiments obtain a variety of stable steady states on larger systems as a function of 
motor density. These include a "lattice of asters" state in which asters are the only stable 
structures, a "lattice of vortices" state, in which individual vortices are stable while asters are 
absent, as well as an intermediate "aster-vortex mixture" state. The LK approach predicts 
that a single vortex should be the stable state at large motor densities even for very large 
systems. Experiments, however, always see a "lattice of asters" in this regime. 

The LK model predicts that motor density profiles in asters are always simple decaying 
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exponentials, independent of the rates at which motors hop on and off the filament. However, 
experiments and theoretical work suggest more complex decays. Such decays include the 
intriguing possibility of power-laws with exponents which vary continuously as a function 
of the on-off hopping rates. Nedelec, Surrey and Maggs (NSM) derive equations for motor 
profiles around a single preformed aster, showing analytically that such profiles are pure 
power-law in nature[19]. The NSM equations describe motors as either freely diffusing in 
the solvent or as bound to and moving along a microtubule. These states are allowed to 
interconvert. The densities of bound and free motors are governed by separate equations. 
LK, in contrast, use a single equation for the full motor density field which is effectively the 
"sum" of these equations but ignore the dynamics of the difference field. 

The NSM equations describe a single aster configuration composed of a fixed number of 
inward (or outward) pointing microtubules. As a consequence, the density of microtubule 
decreases radially outward from the aster core. Since the free motors become bound only 
in the presence of microtubule, the conversion rate should be proportional to the local 
microtubule density. This yields a non-trivial space dependence for this conversion rate, 
which is responsible for the power-law decay of motor densities. However, NSM do not 
address issues of pattern formation. We suggest here, in contrast to NSM, that at the length 
scales appropriate to a coarse-grained hydrodynamic description of pattern formation and 
at the densities of microtubules in the experiments, it is appropriate to ignore fluctuations 
in the local density of microtubule. The relevant fluctuatin g fi eld is then the orientation 



field for the microtubule at fixed density, as in the LK modelj2C 



Detailed information regarding the ordering of microtubules in the presence of molecular 
motors has come from the extensive simulations of Nedelec and collaborators j? , |^ , 3, Q| ■ 
These simulations, performed in a two-dimensional geometry, obtain asters and vortices, in 
addition to relatively disordered configurations in which both asters and vortices are present. 
The best simulations require as many as 19 parameters to be specified; these include fluid 
viscosities, motor diffusion constants, binding strengths and microtubule bending rigidities. 
However, the uncertainties in these parameters are large, often of an order of magnitude or 
more. 

The hydrodynamic approach used here involves far fewer free parameters (22]. Since 
structures obtained in the simulations are appropriate to specific parameter values simulated, 
it is hard a priori, to guess at what patterns would result on changing one or more of them. 
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However, one might expect several microscopic interactions to yield the same hydrodynamic 



description, provided it is general enough 



22| |. Simulating hydrodynamic equations is often 



far faster and more efficient computationally than direct molecular dynamics simulations. 
This makes efficient scans of parameter space possible [2^. Finally, simpler descriptions in 
terms of a small set of model equations enable analytic progress in some limits. 

We summarize our results here: In a regime of parameter space for our model which is 
closest to that for the LK model, we obtain a single vortex as a stable final state for large 
motor densities. Our results here coincide with the LK results. However, in other regimes, 
asters are favoured. A "lattice of asters" state is stabililized in our model through a low- 
order relevant term in the equation of motion for the microtubule orientation. On small 
systems, constraints due to confinement favour a small number of asters, whose number can 
be increased systematically as parameters are varied. We provide qualitative arguments for 
the stability of different patterns in varied regimes of parameter space. 

We calculate the distribution of free and bound motors in asters and vortices obtained in 
our model. Our results for motor profiles about asters differs from both the LK result and 
the NSM one. We derive an exponential decay of bound motor densities away from aster 
cores, modulated by a power-law in which the exponent of the power law depends in a non- 
universal way on dynamical parameters. The associated decay length for the exponential 
can become very large in some regimes of parameter space, yielding what would appear to 
be pure power- law decays close to the aster core. In contrast, LK would predict a pure 
exponential decay, whereas the NSM result is a pure power law. For vortices, we obtain 
results equivalent to the LK results. 

We obtain, numerically, the solutions to our equations when "crowding" effects due to 
the interactions of motors moving on microtubules are accounted for in a simple way. Such 
effects distribute the motor density more uniformly along the microtubules. We argue that 
the inclusion of such "crowding effects" should further act to favour asters over vortices 
in finite systems. We suggest that the term "spiral" as opposed to "aster" or "vortex" is 
a more generic description of the microtubule configuration. We adduce simulational evi- 
dence for such spiral structures favoured by confinement and point out that the microtubule 
configurations seen in experiments do resemble spirals in many cases. 

Our results relating to pattern formation in systems where the effects of boundaries can 
be neglected appear elsewhere. To summarize these briefly, we show in Ref. that our 
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FIG. 1: Qualitative "phase diagram" illustrating how different states, the disordered state, the 
aster-vortex mixture state, the lattice of vortices state, the single vortex, the lattice of asters, 
dominate in different regimes of parameter space; for a definition of parameters see text. The 
parameter e is plotted on the y axis, with the total motor density m plotted on the x axis. The 
parameter S extends out of the plane of the figure. 

model generates all the patterns seen in experiments, such as the aster-vortex mixture, the 



2J]. We show how these states are linked 



"lattice of asters" and the "lattice of vortices" 
in a non-equilibrium "phase diagram", reproduced in Fig. 1, demonstrating how smooth 
trajectories in parameter space can connect the states observed, in agreeement with the 



experiments |25l| . Fig. 1 is relevant to experimental data in that it shows how a relatively 
small number of parameters may suffice to fix the macroscopic state of the system. 

The outline of this paper is the following: In Section II, we describe the details of our 
model. We discuss and justify each term which appears in our equations of motion for micro- 
tubule and motor fields. Section III describes our numerical methods. Section IV presents 
results from our numerical simulation of the equations of motion for different boundary 
conditions. In Section V, we present analytic results for the profile of motor densities in 



single vortex and aster configurations. In a subsection, these analytic results are compared 
with results from direct simulations of the hydrodynamic equations. Section VI discusses 
how saturation effects resulting from the interactions of bound motors affects motor density 
profiles. Section VII combines the results of Sections V and VI in a discussion of the relative 
stability of aster and vortex configurations. We provide a simple intuitive argument for the 
stability of different patterns as parameters in our model are changed. The concluding sec- 
tion, Section VIII, summarizes the results of this study and outlines possibilities for further 
work. 

II. MODEL 

Our model treats motors attached to microtubules differently from motors which diffuse 
freely in solution. Motors which move on microtubules are referred to as "bound" motors, 
while those which diffuse in the ambient solvent are referred to as "free" motors. These 
are described by coarse-grained fields denoted by and m/ respectively and obey different 
equations of motion. In the absence of interconversion terms changing a bound motor to 
a free motor, rrih obeys a continuity equation involving the current of motors transported 
along the microtubules. The free motor field rrif obeys a diffusion equation with a diffusion 
constant D. These two fields are coupled through mechanisms which convert "free" motors 
to "bound" motors and vice versa. The equations they obey have the form 

dtTTif = DV^TTif - if^ipif + 7b^/m6 (1) 
dtrrih = -AV.{miT) l'f^i,mf - ^'b_,fmb (2) 

j'f^b and 7b^/ are the rates at which free motors become bound motors ("on" rate) and 
vice- versa ('off" rate). These rates have dimensions of inverse time. Note that the total 
motor density field m — mb{r) -|- m/(r) is a conserved field. 

The term AV ■ {uibT) describes the motion of bound motors along microtubules with 
velocity A. One could, as we will do later, make the velocity itself density dependent. (This 
is natural when modeling motor motion in terms of the one- dimensional driven diffusive 
dynamics of interacting particles for which the current J(p) is a non-linear function of 
the density p.) While the rate 7/_>6 should depend on the local density of microtubule, 
we will assume that microtubule density fiuctuations are suppressed at the scales relevant 
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to a hydrodynamic description and retain only the orientational degree of freedom of the 
local microtubule field. We work in two dimensions throughout, since the experiments were 
performed in a quasi-two dimensional geometry jl^. 

The dynamics of the microtubules, given by the equation below, incorporates the terms 
used by Lee and Kardar in specific limits. It also includes one completely new term. As in 
the LK model, we ignore fluctuations in the density of microtubules, concentrating on their 
orientational degrees of freedom. The hydrodynamic equation includes terms which reflect 
the dynamics of individual microtubules. We take these to be stabilized at unit length. It 
also includes motor-independent and motor-dependent orientation terms. In principle, for a 
non-equilibrium problem, all symmetry allowed terms must figure here. Of these terms, we 
will incorporate only the lowest-order symmetry allowed terms whose contributions can be 
justified transparently on physical grounds. 

Our equation then reads, 

dtT = T{a - (3T^) + -frribV^T + iVm^ ■ VT + k'V^T + S'Vmb (3) 

The first term, T(a— /?T^), governs the stabilization of the microtubules at a preferred length 
of which we will normalize to unity. The second and third terms, 7m;,V^T+7'Vmfe- VT, 
are alignment terms, reflecting the alignment of microtubules due to the action of bound 
motors. The third term is also interpretable as a "convective" term, in which the local 
velocity which convects fluctuations in the T field is proportional to the gradient of bound 
motor density. 

The fourth term, /t'V^T, describes an intrinsic stiffness against distortions, allowing 
tubules to form an ordered phase in the absence of thermal noise. (This term was dropped 
by LK, where it was assumed that the orientation of the tubules was solely driven by the 
motor current. We include it here because it is certainly symmetry allowed, is linear in the 
field T and thus should appear at lowest order. It is also useful in maintaining the local 
smoothness of the patterns generated). 

The last term is completely new; it is a symmetry allowed term of linear order in the 
fields. Note that a few terms allowed by symmetry have been intentionally excluded. These 
include the term m;,V(V ■ T); such a term has an effect equivalent to the Vm;, term of the 
equation above which is lower order in gradients. 

In the LK model the coefficients 7 and 7' are taken to be equal. Thus, the second 
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and third terms can be interpreted in terms of the functional derivative of a "free energy" 
term. In general however, away from thermal equilibrium, these two coefficients differ and 
we may explore the regimes of parameter space in which their relative strengths vary. For 
convenience we choose 7' = e7 and vary the parameter e to tune the ratio of these two terms. 
The following transformation simplifies the equations considerably: scale length in units 



of D / AJ [3 / a, time in units of [3D / {aA^), motor density in units of D/7 and the tubule 



density in units of \Ja/ [3. The equations then reduce to 

dtruf = V^nif - Jf^bmf + Jb^fnii, (4) 

dtfUb = -V-inibT) + 'jf-.bmf - ^b-^fUib (5) 

dtT = CT(1 - T^) + nibV^T + eVnib ■ VT + kV^T + SVrrib (6) 

The parameter C given by jSD/A"^ is the growth constant. 7/-»f, and 7;,-+/ are scaled in 
units of inverse time (i.e. (/3D/(a;A^))^^. k' is appropriately scaled to k, = n'/D and 
S = S'{l3D)/{a'yA). Note that the scaled equations for the free and bound motors (Eqns. 
I4I5|) are invariant when the motor densities are multiplied by a constant. We will use this 
invariance in Section V(C) to compare analytic results for motor density profiles with results 
from numerical simulations. 

We relate our scaled parameters to typical experimental values in the following way: The 



tubule density, scaled in terms of ya//3, is chosen to be unity. The diffusion constant D 
is about 20/im^/s and A ~ 1/im/s, defining basic units of length and time as 20/xm and 
20 seconds respectively. A tubule density of 1 implies that over a coarse-graining length of 
400/im^, there are around 400 microtubules, a value close to that used in the simulations isl ]. 
Our cho,ce for ov-. a.d conespo^ds to phys.cal rates of 0.005.- to 0.05.-. shghtly 
smaller than those in the simulations |16l|: using larger rates does not affect our conclusions 
here. 



III. NUMERICAL METHODS 

We solve Equations I4I5I and lUl numericallv on an L x L square grid indexed by with 
i = 1, . . . L and j = 1, . . . L. The equation for the free motor density is evolved through an 
Euler scheme, 

mf{t + At) = mf{t) - AtV ■ J{mf) - '^f^bf^if + '^b-^f^b-, (7) 
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where 



Jx{mf{i,j)) = {mf{i + - mf{i - l,j))/2S, (8) 
JyMJ)) = (rrifiij + 1) - mf{i,j - l))/25. (9) 

The grid spacing is 6x = 6y = 6 = 1 and the time step At = 0.1. At the boundaries, we 
impose the boundary condition that no current (either of free or bound motors) flows into 
or out of the system. This condition is easily imposed by setting the appropriate current to 
zero. 

A related discretization is used for the bound motor density equation 

dtmb = -V • (mfeT) + ^f-^trrif - jb-^fmb, (10) 

where, in the bulk, partial derivative terms are discretized as 

d,{mbT) = {mb{i + l,j)T,{z + 1, j) - mb{i - l,j)T,{i - l,j))/25. (11) 

with a similar equation used for the y-component. 

The T equation is differenced through the Alternate Direction Imphcit (ADI) operator 
splitting method in the Crank-Nicholson scheme. At the first half time-step 

(rl - JC,)Tr'/' = (rl + £,)r;, (12) 

where r = 25^/ 6t and I is the identity matrix. The superscripts on Tr^ indicate the time-step 
at which these quantities are calculated. The operators and Cy are given by 

jC^ = rrihdl + {d^mb)d^ + S{d^mb), (13) 
Cy = 2C(1 - T") + rubd^ + {dymb)dy + S{dymb). (14) 

The first and second derivatives evaluated for a function f{i,j) on a lattice point (i,j) in 
the bulk are 

dj{i,j) = {fit + 1, j) - fit - 1), i))/25, (15) 
dlfii,j) = (/(^ + l,J)-2/(^,J) + /(^-l,J))A^ (16) 

with similar equations for the ^/-derivatives. 
At the second half time step 

{rl - Cy)T:+' = {rl + C.)T:+'/', (17) 
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where 



Cy = rribdl + {dymb)dy + S{dymb). 



(18) 
(19) 



A similar scheme is used for differencing the equation for the Ty component. Our simu- 
lations are on lattices of several sizes, ranging from L — 30 to L — 200. We vary the motor 
density in the range 0.01 to 5 in appropriate dimensionless units. We work with two differ- 
ent types of boundary conditions on the T field. In the first, which we refer to as refiecting 
boundary conditions, the microtubule configuration at the boundary sites is fixed to point 
along the inward normal. In the second, which we refer to as parallel boundary conditions, 
microtubule orientations at the boundary are taken to be tangential to the boundary. In 
both these sets of boundary conditions, the state of the boundary T vectors is fixed and does 
not evolve. The total number of motors, initially divided equally between free and bound 
states and distributed randomly among the sites, is explicitly conserved. 

We add weak noise, primarily in the T equation of motion, to ensure that true steady 
states are reached in our simulations. Large noise strengths wipe out patterns, yielding 
homogeneous states. We have also experimented with a variety of initial states to ensure 
that the quahtative features of the patterns we obtain as stable steady states are, indeed, 
robust. 

IV. RESULTS AND DISCUSSION: NUMERICAL 

This section presents the results of our simulations in different regimes of parameter 
space for the boundary conditions discussed above. The simulations discussed here are on 
systems of size L — 30, corresponding to physical length scales of about 60//m. Our results 
on somewhat larger systems (L = 50) are intermediate in character between the results for 
L = 30 presented here and our results on systems of much larger sizes (L > 100), where 
boundary effects are negligible. 

We work at a fixed large motor density of m = mi, + nif = 0.5. We work at fixed values 
of 76-^./ = 7/->6 = 0.5 here but have checked that making the motors more or less processive 
does not alter our results qualitatively. For motors with very small "on" and "off' rates, we 
see disordered states best described as aster-vortex mixtures. The patterns obtained here 
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FIG. 2: Steady- state configurations of microtubules with reflecting boundary conditions. The 

parameters are rn = 0.5 and 7/^6 = 76^/ = 0.5. Patterns are shown at different e at (a) 6' = 0, 
{h)S = 0.05 ,{c)S = 0.5 and (d) S = 2. 

emerge at still higher motor densities for low values of the processivity. We work with a 
small value of k, {k — 0.05 here), such that the magnitude of the self alignment term kV^T 
is small compared to rribV^T. 

A. Reflecting boundary conditions 

Our results for reflecting boundary conditions are shown in Fig. 2(a) - (d). Each column 
represents a different value of S : i.e. (a) S = 0, {h) S = 0.05 (c) S = 0.5 and (d) S = 2. 
We vary the parameter e in these scans, with e taking the values e = 0.0, 0.5, 1.0 and 5 as 
shown. 
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FIG. 3: Steady- state bound motor density profiles with reflecting boundary conditions. The 
parameters are m = 0.5 and 7/^6 = Jb—^f = 0.5. Profiles are shown at different e at (a) S = 0, 
{h)S = 0.05 ,{c)S = 0.5 and (d) S = 2. The darker regions indicate regions of lower motor density. 

In Figs. 2(a) - (d), the effects of varying e with reflecting boundary conditions are the 
following: In Fig. 2(a), S — and the steady state configuration at e = is an aster 
with tubules directed radially inward towards the core. When e is increased to 0.5, the 
configuration resembles a "spiral" . As e is increased further the spiral distorts into a vortex. 
In Fig. 2(b), we show configurations at small but non-zero S, S = 0.05. The single aster is 
stable for e = and 0.5, but yields to a single vortex for larger e. A general observation is 
that the core of the vortex is increasingly distorted as e is increased further. 

Fig. 2(c) shows steady state configurations at 5" = 0.5. A tendency towards the formation 
of a lattice of asters is apparent, as e is increased at these values of S. This is consistent with 
our earlier results on large systems, where we observed that non-zero S always promoted the 
lattice of asters. At large e e.g. e = 5, there is a pronounced tendency towards alignment, 
yielding a steady state pattern of a vortex with a highly distorted core. Fig. 2 (d) shows 
configurations at 5* = 2, where the lattice of asters is present at all values of e shown. At 
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FIG. 4: Steady- state configurations of microtubules with parallel boundary conditions. The 
parameters are m = 0.5 and 7/->6 = 7&->/ = 0.5. Patterns are shown at different e at (a) S = 0, 
(h)S = 0.05 ,{c)S = 0.5 and (d) S = 2. 

large e, the tendency towards parallel alignment competes with the tendency towards aster 
formation. At still larger values of e, we obtain an aligned phase in the bulk (not shown), 
reminiscent of the "bundles" seen in the experiments at large motor concentration. 

Figures 3(a) - (d) show bound motor density profiles corresponding to the microtubule 
arrangements of Figs. 2(a) - (d). Lighter regions of the figure indicate regions of larger 
density. Note that the bound motor density is concentrated at the centres of asters. In 
vortex configurations, the motor density profiles for bound motors are far smoother. In the 
"lattice of asters" configuration, the bound motor density profiles peaks at the centres of 
the asters, decaying to a small value at intermediate points far away from aster cores. The 
profiles of free motor densities are visually very similar to those for bound motors and are 
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FIG. 5: Steady- state bound motor density profiles with parallel boundary conditions. The 

parameters arc m = 0.5 and = 76-+/ = 0.5. Profiles are shown at different e at (a) S = 0, 

{h)S = 0.05 ,(c)S' = 0.5 and (d) S = 2. The darker regions indicate regions of lower motor density. 

not shown here. 

B. Parallel boundary conditions 

We now discuss pattern formation in finite systems with parallel boundary conditions, see 
Figs. 4(a)-(d). These are to be contrasted to Figs. 2(a) -(d). Note that qualitatively differ- 
ent sequences of patterns are stabilized at identical values of other parameters, depending 
on the boundary conditions. This illustrates the sensitivity to boundary conditions which 
obtains for sufficiently small system sizes L. 

Fig. 4(a) shows pattern formation with parallel boundary conditions at -S" = 0. A single 
vortex is obtained as a steady state at all values of e shown i.e. e — 0.0, 0.5, 1.0 and 5. Fig. 
4(b) shows patterns for S — 0.05. The steady state here is a single aster with tubules pointing 
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outwards at the core, but aligned with the configurations at the boundary. (We term such 
configurations as "outward asters".) Clearly these arise as a consequence of the boundary 
conditions which favour vortices and the parameter regime which favours asters. At e = 5, 
the steady state is a well-formed, clean vortex. Fig 4(c) illustrates pattern formation at 
S = 0.5, as a function of e. We see a tendency towards the formation of the lattice of asters 
phase expected for large S for sufficiently large e, although indications of this are seen even 
for small e. 

At very high e, we obtain an unusual configuration which we term a "fiag". In this 
configuration, tubules point radially outward near the core while merging with the imposed 
parallel microtubule configuration near the boundary. However, the configuration appears 
to have a four-fold axis, in which the axis lines are along the diagonal. Across this axis, 
the microtubule orientation changes sharply. This 4-fold symmetry refiects the four-fold 
symmetry of the simulation box. Finally, Fig. 4(d) shows patterns formed at 5" = 2. In this 
regime, the lattice of asters is the stable steady state for all values of e. 

Figures 5(a) - (d) shows bound motor density profiles corresponding to the microtubule 
arrangements of Figs 4(a) - (d). Note that the density distribution in single vortices varies 
smoothly, consistent with theoretical expectations. The profiles appear sensitive to the 
boundary and a four-fold rotation axis can be seen in several of the patterns which involve 
a single vortex. The patterns for non-zero and large ^S" are very similar to those obtained for 
refiecting boundary conditions with similar values of parameters. The profiles of free motor 
densities are again visually very similar to those for bound motors and are not shown here. 

C. Discussion 

Our results for the cases outlined in the previous sub-section are summarized as fol- 
lows. For refiecting boundary conditions, we obtain the general sequence aster — > spiral — > 
vortex at yS" = 0. This sequence is obtained at fixed motor density, as a function of the 
non-equilibrium parameter e. With parallel boundary conditions, the patterns formed are 
generically vortices, although the region surrounding the core is progressively distorted as e 
is increased. Finite S favours the formation of a lattice of asters, as in large systems. We 
observe some novel configurations such as the "flag" configuration, refiecting the four-fold 
symmetry of the simulation box and the outward aster. 
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We commented earlier that the configurations obtained in experiments appear to repre- 
sent states intermediate between asters and vortices in the most general case. We expand 
on this here: A spiral configuration is described by the unit vector field T = T^f + Tg9 in 
which Tj. and Tg take the forms, 

Tr = cos(a) , Tg = sin(a), (20) 

where a is a constant. This contains both asters and vortices in appropriate limits: In the 
limit a = 0, this equation describes an aster, while in the limit a = 7r/2, the configuration is 
a vortex. Thus asters and vortices are particular limiting cases of more general spiral states. 

In our simulations, spirals arise as a consequence of the competition between the vortex 
configurations favoured locally for intermediate and large values of e and refiecting boundary 
conditions, which favour asters. Spirals offer the best compromise between these and it is 
reasonable to expect that such states should be more generically seen than either asters or 
vortex states. Indeed, inspection of the configurations associated to the "lattice of vortices" 

n 

in Ref. jl3|] provides strong visual evidence for generic spiral states. 

The bound and free motor profiles we obtain in the simulations are consistent with 
expectations from our theoretical analysis (see below). The characteristic sharp peaks in 
motor density profiles associated with asters are replaced by far more slowly varying profiles 
for vortices. The lattice of asters, therefore has strong signals in the motor distribution 
function. Appropriate experimental labelling of motors in this phase should yield patterns 
and profiles closely similar to those displayed here. 



V. ANALYTIC RESULTS FOR MOTOR PROFILES IN ASTER AND VORTEX 
CONFIGURATIONS 

We now present an analytic calculation of free and bound motor density profiles in vortex 
and aster configurations. We fix single vortex and aster configurations of microtubules and 
solve for steady state free and bound motor densities in the presence of the fixed microtubule 
configuration. 

We find steady state solutions of Eqns. El and El by setting the time derivatives to zero 
i.e. dtTTif = dtrrib = 0. Adding and subtracting these equations, we obtain 

V^nif - V • (mfeT) = 0, (21) 
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V'^rrif + {'^b^fmb - 'jf^brrif) = 0. (22) 
The first of tliese equations implies tliat 

Vm/ = rribT + C, (23) 

wfiere C is an arbitrary vector wliose divergence V ■ C = 0. 

A. Single vortex i.e. T = 9 

Assuming radial symmetry, we may write mj = mf{r) and nib = mb{r). Let C = Crr+ceO. 
Eqn. 1^ then implies 

drtn ff = nibO + CrT + cqO . (24) 

Hence, 

Cr = drTTif and C0 = —nib. (25) 
The constraint of zero divergence for C yields 

drCr + ^ + ^deC0 = 0. (26) 

Since cg = —mb{r) is purely radial 

dece = 0. (27) 

Eqn. |2ni then becomes a radial equation 

drCr + - = 0. (28) 
r 

Using Cr = drTTif, we obtain 

V^m/ = 0, (29) 

which supports solutions of the form 

m/(r) = ci + C2 ln(r), (30) 

where Ci and C2 are constants to be determined by boundary conditions and normalization. 
Eqn. 1221 now simplifies to 

jb^frub - If^brnj = (31) 
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which yields 

mbir) = ^{ci + C2Hr)) (32) 

We must determine the constants Ci and C2. There is no radial current of motors (bound or 
free) in the vortex configuration. The condition of no radial current of free motors at the 
boundary implies 

drmf{r) — ai r — L, (33) 
therefore implying that C2 = 0. The constant Ci is determined by normalization 

J £r(mf(r) + mb(r)) = N. (34) 



In the vortex configuration 
which gives 



/ ciV(l + ^)ci = TV, (35) 



^1 = + (36) 



with n the density of motor N/ L'^. 

The result above holds for purely tangential boundary conditions, consistent with the 
symmetry of the vortex. It will be modified vis. a vis the simulation results both as a 
consequence of the (cartesian) latticization used to discretize the equations as well as by the 
four-fold symmetry of the simulation box. Both these factors will lead to a combination of 
the logarithmic and constant solutions above. 



B. Single aster solution i.e. T = — f 

We again assume radial symmetry for the bound and free motor densities. We choose 
C = f{r)r, with /(r) such that C is divergenceless i. e. V • C = 0. We thus obtain 
/(r) = Co/r, with cq a constant. The radial current of free motors is —drUifr while the 
radial current of bound motors is —mhf. The boundary condition that the total motor 
current vanishes at the boundary implies 

drTUf + mb — co/r = at r — L. (37) 

This ensures that Cq = and therefore 

drmf{r) = —mb{r). (38) 
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Eqn. 1221 becomes 



d'^nif + - jb->f)drmf - jf^brrif = 0. (39) 



This is a second order differential equation whose solution is completely specified if the value 
of the function and of its derivative at the boundary are supplied. Given these boundary 
conditions, the free motor density can be obtained via a numerical solution using the Runge- 
Kutta method. The bound motor density is also easily determined via a radial derivative of 
the free motor density. 

We now derive exact and asymptotic expressions for motor densities in the aster geometry. 
The general solution to the equation above is a combination of confluent hypergeometric 
functions and has the form 

. 1 ib^f 



, ,h^lLf + ^V^i.r)]. (40) 

2 2J^f!Z;+4^b 



The functions ,F,{1 - .^^^^^^^^^^^ , 1, V^tL/ + 47/^6 r) and f/(l - 
1, A/Tb^/ + 47/^6 r) are the two solutions of the confluent hypergeometric 



7i>->/ 



Kummer equation. It is useful to define a quantity p given by 

p = ki- )■ (41) 

Note that 

< p < 0.5, (42) 

with 7b^j,7j^;, > 0. The argument of the exponent in Eq. EUlis always negative for all 
7b^/,7/^fe > 0. The coefficients Ci and C2 are determined by boundary and normalization 
conditions. 

The boundary condition is that there is no total motor current at the boundary, i.e. 

drrrifir) + m{,(r) = at r = L, (43) 

since the radial current of bound motors in the aster configuration is J = —mh{r)f. This can 
be used to determine C2 in terms of Ci. We find that on imposing the no-current boundary 
condition, C2 is very small compared to Ci and is significant only when the system size is of 
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order the "correlation" length over which the density decays. We will assume (see below) 
that we can set C2 = to yield physically admissible solutions. 

The constant Ci can now be fixed by the normalization condition which requires the total 
number of motors (bound+free) be constant, 

/ d\{mb{r) +mf{r)) ^ N. (44) 

t/ 

where is the total number of motors. 

A physical argument for neglecting the iFi{^ — — ^=2^s4^=, 1, ^J'y^^f + 47/^6 r) term 
in comparison with the — — , 1, sj'yb-^f + ^lf-*h f) is the following. A power 

series expansion yields 

where (p)„ = (p + n — l)!/(p — 1)! and (p)o = 1. This function is a monoton- 
ically (exponentially) increasing function of r. Such a density distribution implies 
a motor density per unit area which increases as r is made arbitrarily large, a so- 
lution which is clearly untenable on physical grounds. On the other hand, U{\ — 



^-^^5^^==, 1, y^7^_^j + 47/^6 r) decreases monotonically with increasing r. Hence we 

retain only the U{\ / ^"'^^ 1; Jlh-^f + 47/^6 r) part of the solution for mf{r). 

To derive the asymptotics we begin with the integral representation 

U{p, Q,^) = Y^^ I e-^H^-\l + ty-^-'dt. (46) 

We relabel + 47/_>b)r — z and |(1 — ^ 2 ^''^4 =) — P notational convenience. 

We then obtain 

1 roc p-zt j.p-1 

U(p, 1, z) = —— / —dt, (47) 

We now change variables to u = zt. The integral is then 

1 1 .00 e-"MP-^ 



^/ X / / , , du. (48) 
r{p)zPJo {l + u/z)P ^ ' 

Since p is a small number between and 0.5 we may expand the denominator binomially. 
The first term gives r(p) and subsequent terms converge. Hence in the large z (asymptotic 
hmit), the integral is just r(p). Therefore, 



mf{r) = ci 
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FIG. 6: Density profiles for bound motor densities. The profiles are shown for an aster in the 
regime 7/-+6 = 76->/ = 0.5, with S = 0. The density is chosen to be 0.5. 

where the inverse of the "correlation" length is defined as 



The correlation length ^ and the power- law exponent p depend on 7/_>b and 76_»/. We see 
that the bound motor density in the aster case has an exponential fall modulated by a 
power-law tail instead of the pure exponential falloff predicted in the LK model. 

Another limit in which an exact answer can be obtained is the following: If there is no 
interconversion between the two species of motors {i.e. 'Jb^f = If^b = 0), the equations of 
motion for free motors and bound motors decouple. The number of free motors and bound 
motors in the system are then conserved independently. The free motor density equation 
supports solutions of the form m/(r) — ci + C2ln(r), where ci and C2 are constants of 
integration. Imposing a vanishing current of free motors at the boundary constrains C2 — 0. 
The other constant Ci can be fixed from the condition that the total number of free motors 
is constant. The bound motor density then has the simple power-law behavior mb{r) ~ 1/r. 




(49) 
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FIG. 7: Density profiles for free motor densities. The profiles are shown for an aster in the regime 
7/_+f, = 7b-»/ = 0.5 and 5 = 0. The density is chosen to be 0.5. 

C. Comparison of the analytic results with simulations 

These analytic results can now be compared to the results of direct simulations of the 
dynamical equations for single aster and vortex structures. We choose T = f in the full 
equations of motion for the motor densities given in Eqns. ^ and (0) to represent the aster. 
We then evolve these equations in time till the steady state is reached. We compare the 
motor profiles obtained in such simulations to profiles obtained from solving the ODE at 
steady state (Eq. EH)) by a 4*^ order Runge-Kutta method. Free and bound motor densities 
at the boundary of the system are obtained by solving the full time-dependent equations 
(Eqns. Eland EI). These are used as boundary value input to the Runge-Kutta procedure to 
facilitate direct comparison between simulation and analytic results. 

Using the symmetry of the scaled equations Eqns. |3] and El discussed in Section II, we 
normalize the densities obtained by these two methods to the same constant value before 
comparing them. Figure 6 shows a plot of the distribution of bound motors in a single aster 
obtained from (i) the direct simulation (ii) the Runge-Kutta method and (iii) the exact 
solution outlined in the previous section. The parameter values are 7/^6 = 76^/ = 0.5 and 
m = 0.5. Fig. 7 shows the comparison between simulations, the Runge-Kutta method and 
the exact solution for the free motor density profiles. While the data is noisy, particularly for 
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FIG. 8: Density profiles for bound motor densities. The profiles are shown for a vortex in the 
regime 7/^6 = Ib^f = 0.5 and 5 = 0. The density is chosen to be 0.5. 

the bound motor densities, the agreement with the theoretically predicted result is satisfying. 

Fig. 8 shows the distribution of bound motors in a vortex configuration obtained by 
fixing T = in the full equations of motion for the motor densities given in Eqns. lU and El 
The plot, of mf,{r) vs. r, illustrates the slow, essentially logarithmic variation of the motor 
densities in such configurations, consistent with the analytic approach of this paper and of 
earlier work. 

VI. SATURATION EFFECTS AND THE STABILIZATION OF ASTERS 

So far, in our analysis, we have ignored the effects of interactions between bound motors. 
These motors, in their physical setting, have a finite size and move on a one-dimensional 
track, the microtubule. It is thus reasonable to expect that interactions between bound 
motors should become increasingly important at large motor densities. One can account for 
these interactions by simply requiring that more than one motor cannot occupy the same 
volume in space at the same time. Such "excluded volume" interactions have been used in 
earlier models for collective effects in molecular motor transport [27 1 . 

At the simplest level, accounting for such interactions leads to a non-linear (in general, 
also non-monotonic) dependence of the motor current on the density. This can be incorpo- 



• Simulations 
5.3 - 2log(r) 
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rated in our calculation by choosing a motor current of the form 



= Am,f{mh)T. (50) 

The function f{mi,) should, in principle, be calculated from an underlying microscopic model. 
We will circumvent this necessity by simply assuming a convenient form purely for illustra- 
tive purposes which is (a) consistent with the requirement that A (the velocity) is density 
independent for small rrih densities and (b) saturates for large rrib. One such choice is 

f{mb) = 1 - tanh(mb/msat). (51) 

The value rrisat limits the current of bound motors. When rrisat ^ f^b then f{mh) — > 1 and 
we recover results discussed earlier. When niaat ^ "^6, the current reduces as f{mi,) 0. 
In the aster case (T = -f), Eq. IHHl becomes 

drmf{r) = -mb{r) f {mb) . (52) 

Therefore, the free motors obey 

d^rrif + - 1b^f/f{mb))drmf - Jf^brrif = 0. (53) 

This equation can now be solved as a boundary value problem using the methods described 
earlier. 

Figure 9 shows bound motor densities mb{r) as a function of r for different values of 
'nT'sat- We have normalized each of the plots to the same value at the origin. We choose 
iiisat = 0.001,0.1,1 and 1000. An increased rrisat implies smaller saturation effects. The 
plots shown in Fig. 9 are normalized to unity at r = 1. Observe that reducing nisat leads 
to far smoother variations of the bound motor density. Data for the associated free motor 
densities are similar and are not shown here. 

VII. RELATIVE STABILITY OF ASTERS AND VORTICES 

As we have seen in Section IV, single asters are unstable to single vortices for large 
systems, provided S = 0. This is surprising because both asters and vortices are equivalent 
topological defects - an aster can be transformed into a vortex by rotating all T vectors 
through 90°. However, motor density profiles in aster and vortex configurations are different. 
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FIG. 9: Effect of saturation on the motor density profile. As the effect of saturation is increased, 
the profile is smoother. The curves are labeled from left to right as mgat = 1000, 1.0,0.1,0.001. 

It is thus natural to seek an explanation for the relative stability of asters and vortices at 
S — which proceeds from this observation. 

Lee and Kardar rationalize the stability of vortices vis. a vis asters in the simulations in 
the following way: If e = 1, the right-hand-side of the tubule equation derived by LK can be 
interpreted as the functional derivative of a "free energy". The motor density profile plays 
the role of a spatially modulated stiffness. This enables the calculation of the relative "free 
energies" of aster and vortex states. Smoother motor density profiles lead to lower "free 
energies" . The reduced energy of the vortex configuration, in which motor density profiles 
decay logarithmically, implies its increased stability as compared to asters, for which motor 
density profiles decay exponentially away from the core. 

The competition between strain energies of asters and vortices thus implies a system-size 
dependence of "energies" and a cross-over length scale at which the energies of the two are 
comparable. For system sizes smaller than this length, asters are favoured while vortices are 
favoured for larger system sizes. This simple argument rationalizes the relative stability of 
aster and vortex states. We will use it to explain an intriguing feature of our data - the 
presence of stable asters for fairly large system sizes e.g. L = 50, in the S = limit. At 
these scales and for equivalent parameter values, the LK model would have yielded a single 
vortex in steady state. 
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FIG. 10: A plot of the difference in energies of a single vortex and a single aster versus the system 
size L. The parameter values are m = l,e = 1 and S = 0. The y-axis is plotted with a logarithmic 
scale. The point at which the energy of the vortex becomes lower than the energy of the aster 
shows up as the minimum in this curve. In the Lee-Kardar model the length scale at which the 
vortex becomes lower in energy is about four times smaller than in our model. 

We adapt the Lee-Kardar argument to our model, using our analytic results for bound 
motor density profiles in asters and vortices. In Fig. 10, we show the modulus of the difference 
in energy of a single vortex and a single aster as a function of the system size for our model 
and for the Lee-Kardar model. We work at the same value of total motor density m in 
each case, setting e — 1. Since the magnitude of "energy" scales in the two models differs, 
the behaviour of the energy difference is most clearly seen using a logarithmic scale for the 
y— axis. The point at which the vortex energy is reduced below the aster energy appears as 
the minimum in the curve. From the position of this minimum, we see that the crossover 
length scale above which asters are unstable to vortices is smaller in the Lee-Kardar model 
as compared to our model. 

The fact that this crossover lengthscale hes between 50 and 100 in our units for length 
rationalizes our observation that the stable asters we obtain for L — 50 are replaced by stable 
vortices for L = 100. These results support, and add further credence to, the Lee-Kardar 
argument, at least in the 5 = case. With finite S, of course, vortices are disfavoured 
altogether and generating stable aster structures is no longer an issue. 
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We now provide a qualitative explanation of three basic features of motor and microtubule 
arrangements at zero and non-zero S. These are (i) the absence of vortices at non-zero S, 
(ii) the presence of a small number of large asters at small S vs. (iii) the presence of a large 
number of small asters which increases as S for larger S. Our approach is approximate and 
qualitative, based on "free energy" type arguments, analogous to those of Lee and Kardar. 

Consider the limit in which e — 1 and the Vm term is obtained through a functional 
differentiation of the "free-energy" term 

Es= [ (fxiaS/ ■ T + -mb)^ (54) 
J a 

A functional derivative of F with respect to the T field generates the terms a^V(V ■ T) and 
SVmi,. In the limit that a — >■ 0, we obtain the Vmj term of our equation of motion. (We 
prefer to work with a non- vanishing, although infinitesimal, value of a and to consider the 
effects of varying to avoid the singular behaviour which arises in the a — Q limit.) Our 
introduction of the term in S was, in fact, motivated by the observation that motors moving 
along two initially parallel microtubule configurations, act to bring them together, favouring 
a non-zero divergence. This physics, though incorporated in the simulations of Nedelec et 
ai, is ignored in the LK model. 

The full "free energy" for the T field is 



E= (fx 



■la\T\' + ]f3\T\' + lm,{r){VTf + haV ■ T + -m,f 
z 4 z z a . 



(55) 



We fix the modulus of T, thereby eliminating the first two terms, and concentrate on the 
relative energetics of states governed by the last two terms of the equation. 

For the Lee-Kardar model, with S = 0, the stability of vortices over asters can be 
understood from simple comparisons of the relative energies of vortices and asters. First 
consider the case of a single aster. We choose a simplified profile for bound motors, assuming, 

rrib — A r < ^ 

mfe = r>^. (56) 

This represents, although approximately, the fact that microtubule densities are significant 
only over distances of order ^ from the core of the aster. The factor ^4 is a normalization 
factor which ensures a fixed number of motors in a system of typical size L. We will work 
with a fixed overall density of motors ~ N/L"^ and examine the limit in which L ^ oo. 
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In the S — limit, the energy Ea of the aster is dictated by 

Ea^ /^ciV m6(r)(VT)2~^log(e/a) (57) 

Ja 

where we cut the integral off at a lower limit a corresponding to a microscopic coarse- graining 
scale and have used the fact that (VT)^ — for both asters and vortices. We now fix 
A by normalization: since m{r)£r ^ N ^ nL'^ we obtain ~ nL^ log(^/a)/^2_ 

(Accounting for a more complex decay of the bound motor density does not change this 
result qualitatively.) Thus, for a system of size L (L ^ ^), imposing a fixed density of motor 
yields a quadratic increase of the energy with system size L. 

We now repeat the same argument with a single vortex. Here, the relatively slower 
variation mh{r-) ~ Ci + C2log(r/A) yields E^ ~ {mb{r) / r'^)rdr ~ log(L/a), upto further 
logarithmic corrections. The prefactor is then independent of system size, leading to an 
overall logarithmic increase of energy with system size, a far weaker function than the 
quadratic dependence on L of the aster. This imphes that for sufficiently large L, the 
energy of the isolated vortex falls below that of a single aster, as first suggested by Lee and 
Kardar. in systems of sufficiently large sizes. 

How are these arguments modified at finite 5*? Since simulations indicate that the state 
obtained at finite positive 5" is a "lattice of asters" with the scale of the aster decreasing as 
S is increased, we will compare the energies of single vortex configurations with energies for 
an assembly of asters ("mini- asters") of typical size a. We imagine that the Lx L system is 
subdivided into a x a units and place an aster in each one of these sub-units. For a system 
of linear size L, the number density of such asters is . For inward pointing asters, as 
obtained in our calculation, the divergence V • T = — 1/r. Consider first the term 

/• r 5 1^ 

Eg d^x aV • T + -mb(r) . (58) 
J I a i 

Now mini-asters are (i) assumed small and (ii) obey boundary conditions which differ from 
the case of the single system-size spanning aster. Bound motor profiles are then, in general, 
combinations of the two solutions obtained earlier rather than the single one which yielded 
the exponentially damped form used in our earlier analysis. Also, the close proximity of the 
many mini-asters formed indicates a slower than exponential variation of the bound motor 
densities about each one. Let us assume that motor densities in such mini-asters adjust 
in order that the following motor density profile is obtained: rribir) ~ a^/Sr{r < a) and 
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777.6 (r) ~ 0(r > a). This ensures that the contribution of the term above is cancelled. Note 
that the bulk of the contribution to the energy of a mini-aster comes from the region r < ^; 
the contributions from larger regions is negligible provided ^ ~ u. 

We first calculate the energy Ei^n^ of a single such mini-aster configuration. This is 
obtained from the integral 

E,mA = I d'rm,ir)iWTr ^ dr- ^ ^ (59) 

The number of motors in a single mini-aster is obtained from 

NimA ~ / m(r)dV - — ^ (60) 

J a O 



The total energy E'^j^ of the system of miniasters is then 



EIa - ^/L- (61) 



The constraint n — (pNijnA yields 



(62) 



illustrating how the number of asters increases as is increased. The total energy corre- 
sponding to a fixed areal density of motors 77 is then easily derived as 

aa 

We now compare this with the energy of a single vortex at non-zero S, computed assuming 
that the motor density profiles are the same as in the case in which S was zero. This is 
obtained as E^ ~ S'^n^LF' /a^, yielding the ratio 

(64) 



Ev nS'^aa 



We may also compare the energy of an assembly of mini-asters with the energy of a single 

aster. Arguments similar to those above yield Ea ~ t?S'^L^ jo? ^ and thus 

Ea nS'^aa 

w ^ ^^^^ 

Note that both these arguments indicate that as a is reduced towards zero (or S is 
increased from a non-zero vahic), both single asters and single vortices arc unstable to an 
array of mini-asters. At infinitesimal a, provided S is non-zero, our results indicate that the 
number density of mini- asters increases with S (linearly in the simple argument given here), 
with the size of each aster decreasing in proportion. The arguments given here rest on the 
assumption that the motor distribution in asters adjusts so as to minimize the cost for the 
term involving S. No such adjustment can lower the energy of single vortex configurations. 
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VIII. CONCLUSIONS 



This paper presents a hydrodynamic theory of pattern formation in motor-microtubule 
mixtures, accounting for the effects of confinement. We show that the influence of the 
boundaries on pattern formation can be considerable, by illustrating how either asters or 
vortices can be formed depending on how the orientation of microtubules is fixed at the 
boundary. We have explored the parameter space of e, m and S systematically, describing 
the variety of configurations obtained. We obtain density distributions of free and bound 
motors corresponding to the final microtubule configurations. Such plots may be compared 
directly to experiments. We have also compared analytic predictions for motor density 
profiles in isolated vortices and asters with simulation data. Our results presented here 
complement our earlier work on pattern formation on much larger systems, in which the 



effects of boundaries were minimal 



One technical improvement would be to account for variations in the local density of 
microtubule, as opposed to only their orientation. We could then account for the density 
dependence of quantities such as 7/^6. More information from experiments, performed in 
confined geometries using motors with a range of different processivities would also be useful 
in clarifying some of the issues which relate to the simulations described here. It would also 
be interesting to search for the novel configurations we obtain here in different regimes of 
parameter space, such as the "flag" , the distorted vortex and the "outward aster" . 

As we have emphasized here and in earlier work the set of hydrodynamic equations 

we motivate and use allow for a minimal yet complete description of the patterns formed in 
mixtures of motors and microtubules. We have suggested elsewhere that it may be useful 
to think of spindle formation itself as a pattern formation problem which can be modelled 
using continuum hydrodynamical equations in a small number of variables. Further work 
relating to this program is in progress. 

We thank Y. Hatwalne, Madan Rao and G. Date for useful and illuminating discussions. 
PBS thanks IBM for providing computational facilities under a shared university research 
grant. 
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